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ABSTRACT 

Early-type galaxies are characterized by many scaling relations. One of them, the so-called 
fundamental plane is a relatively tight correlation between three variables, and has resisted 
a clear physical understanding despite many years of intensive research. Here, we show that 
the correlation between the three variables of the fundamental plane can be the artifact of the 
effect of another parameter influencing all, so that the fundamental plane may be understood 
as a confounding correlation. Indeed, the complexity of the physics of galaxies and of their 
evolution suggests that the main confounding parameter must be related to the level of di- 
versification reached by the galaxies. Consequently, many scaling relations for galaxies are 
probably evolutionary correlations. 

Key words: galaxies: fundamental parameters - methods: statistical - galaxies: evolution - 
galaxies: formation - 
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Galaxies are huge systems of stars, gas and dust, that have as- 
sembled through a complex history of collapses, internal pertur- 
bations, interactions, mergers and even ejections of some mate- 
rial. This complexity renders quite remarkable the scaling relations 
found between many observables and properties, like size, mass, 
surface brightness, velocity dispersion, magnitude, color, metallic- 
ity, black hole mass or spectral features. The fundamental plane is 
a famous correlation between three variables, the effective radius, 
the central veloci ty dispersion and the s urface brightness within the 
effective radius jDressler et al.l 1 19871 : iDiorgovski & Davisl 1 19871 : 
iNieto et al] Il990h . For more than 20 years of i nvestigations, it 



has resisted a satisfactory phy s ical und erstanding ( Robertson et al 



2006; D'Onofrioet; 



ctory phy: 
afl l2008l: 



Gargiulo et al. 2009; Graves et al 



2009 ; iNigoche-Netro et al.ll2009l : lFraix-Burnet et alj201Ch . 



For (astro)physicists, such scaling relations reveal some un- 
derlying physical laws that are described by a set of equations. 
However, for statisticians, correlations are not always causal, as 
lurking or confounding parameters can do the trick. In general, 
these correlations are called "spurious" although this more specifi- 
cally concerns acknowledged false correlations. 

A famous example of a "spurious" correlation in statistics is 
the relation between ice cream sales and number of drownings. The 
confounding parameter is the high temperature which favours both 
the ice cream consumption and the number of people swimming, 
hence the probability of drownings. So even unrelated phenomena 
can appear correlated. Other examples, like the relation between 
the size of a children and his school level or between the distance 
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of the Pioneer probe to the Earth and the average Earth temperature 
since its launch, have time or evolution as confounding factor. 

Indeed the correlation between children's size and school level 
is statistical: the taller ones being more probably found at higher 
levels. This is reminiscent of a classical reasoning used in the de- 
termination of the cosmic distance ladder: for example, more lumi- 
nous stars are more probably found in bigger galaxies. Simillarly, 
because of evolution, bigger galaxies are more probably more mas- 
sive, more luminous, more metallic, with a higher velocity disper- 
sion, etc. Hence statistically, a lot of properties must appear corre- 
lated, the confounding parameter being the level of diversification. 

In this paper, we show that correlations between the three vari- 
ables of the fundamental plane can easily be the artifact of the effect 
of another parameter influencing all. In other words, we show that 
the fundamental plane may be understood as a confounding correla- 
tion. We propose that the confounding parameter(s) is (are) related 
to evolution. Consequently, many scaling relations of galaxies are 
probably evolutionary correlations. 

In Sect. [2] we derive the conditions for parametric relations to 
yield a planar correlation in a 3-variable space. We then dedicate a 
full section to the fundamental plane (Sect. [3} to show how easily 
these conditions can be fulfilled. We finally discuss why and how 
the evolution of galaxies is probably the confounding parameter 
(Sect.©. 



2 CONFOUNDING CORRELATIONS 

Let us consider three variables Cl t , Cl 2 , ^3 that are all functions of a 
same generic parameter X. For sake of simplicity in this paper, let 
us assume that these functions are power laws: 
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^3 



A\X>' 
A 2 X S 
A 3 X' 



(1) 



with A\, A 2 and A3 being constant. Any linear correlation of 
the form 

log Q. 1 = a log Q. 2 + b log Q. 3 + c (2) 
where a, b and c are constant, is thus expressed as 

p\ogX + logAi = as\ogX + bt logX + alogA 2 + MogA 3 + c. (3) 
This expression should be valid for all X, implying that 



p = as + bt 
logA, = a logA 2 + b log A 3 + c. 



(4) 



This set of two linear equations generally yields solutions for 
a and b. The result is thus a plane in the 3-D space (fli, fl 2 , ^3)- 

The generic parameter X can be a multivariate component, 
making the parametric dependence of Q t , Q 2 and Q 3 mutlivariate. 
For instance, consider power laws with two parameters, Xi and X 2 : 



Cl 2 

n 3 



AiXfXf 
AjXfX*' 

a 3 x;< 



(5) 



Then it is easy to show that a correlation like equation l[2j 
holds if: 



P 
P' 
logAi 



as + bt 
as' + bt' 

a log A 2 + b log A 3 + c. 



(6) 



This set of equations is more constraining than equation ©. 
But if for instance the second relation in equation l[6]l is not exactly 
fulfilled, then the parameter X 2 can be considered as noise adding 
a dispersion to the correlation defined by the two other relations. In 
other words, X 2 could generate a thickness to the plane defined by 
equation (0. 

To be complete, we must discuss the physics in the param- 
eter X. This parameter is supposed to influence all the three vari- 
ables. We distinguish two possbilities: causality or evolution. In the 
first case, the relation between the variables and X is driven by di- 
rect causality, that is some parameter influences directly each of the 
three variables 0.\, £l 2 , £1 3 through physical laws. This is illustrated 
in Sect. 13.21 In the second case, the relation is statistical in the sense 
that each variable is bound to evolve so that even totally unrelated 
variables can show an apparent correlation. This is discussed in 
Sect.[33]and Sect.g] 



3 THE FUNDAMENTAL PLANE OF EARLY-TYPE 
GALAXIES 

3.1 General constraint 

Let us put CI] = r e the effective radius, Cl 2 = cr the velocity disper- 
sion, and £1 3 = L the luminosity, equation {TJ becomes: 



r e = AiXP 
l cr = A 2 X S 
[ L = A 3 X 

The surface brightness p e can be expressed as 

p e = -2.5 log (L/4nrfj + m 

= (-2.5? + 5p)logX + 2.51og(47r) + m 



(7) 



(8) 



(10) 



(11) 



where m is a constant of normalisation. Any linear correlation 
of the form 

log r e = a log cr + bp e + c (9) 
translates to 

plogX + logA[ = aslogX + b(-2.5t + 5p)\ogX + 
a log A 2 + b (2.5 log (47tA, /A 3 ) + m) + c. 

This implies 

p = sa + (-2.5t + 5p)b 
logA { = alogA 2 + Z>(2.5 log (4-7rAj/A 3 ) + rnj + c. 

If a solution can be found for a and b from equation jilt , 
then the equation of the fundamental plane equation {9]l is ob- 
tained. There is no need of any further assumption to explain the 
fundamental plane. This demonstration is made here with a simple 
power-law assumption in equation Q. But this result is also true 
for more complex functions, with equation dl It being replaced by 
more complicated conditions. 

To understand the origin equation (9}, we need to solve equa- 
tion dl It . Since there are too many unknowns, additional input is 
required. Two approaches are possible: either input some a priori 
knowledge to determine the functions of X and derive coefficients 
a and b (Sect l3.2l >. or conversely use the observations to determine 
a and b and derive constraints on the functions of X (Sect |3.3t . 

3.2 A priori physical input 

In this section, we consider the input of a priori knowledge believed 
to be relevant for the physics of galaxies. At first glance, it is quite 
logical to consider that mass, either dynamical, true or stellar, is 
somehow related to the radius, the velocity dispersion and the sur- 
face brightness, essentially because it influences the density of stars 
and their kinematics. So we consider X = M in equation ifTJ. 

One obvious way to link the kinematics to the mass is through 
the virial theorem. Hence, let us assume that the virial equation 
holds and that the ratio between the average squared velocity and 
cr 2 , and the ratio between r e and the gravitational radius, are con- 
stant. This gives 



<r 2 r e oc M => 2s + p = 1 

Using equation i ll It and equation Q2) we obtain 

1 - 2s = sa + (-2.5* + 5(1- 2s)) b 
s(a + 2- 10b) = 25tb -5^+1 

2.5^-56 + 1=0 if a + 2- 106 = 

s = (2.5tb -5b +1)1 (a + 2- \Qb) otherwise. 



(12) 



(13) 



Unfortunately, the brightness has no direct relation to mass, 
but it might be assumed that the mass is essentially due to the stars 
that are responsible for the luminosity. If we also assume that the 
ratio M/L is constant for a given population of galaxies, we must 
have: 



L oc M 



1. 



(14) 



Replacing this value of t in equation Q3), one obtains two 
solutions: 



(0 



2.5b - 5b + 1 = and a + 2 - 10b = 
=> b = 0.4 and a = 2 



(15) 
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■ (2.5b -5b + l)/(a + 2- 10/?) 
*s= (1 -2.5b) I (a + 2- Wb) 



The first solution is well known, being the so-called virial 
plane a = 2 and b = 0.4, a classical solution obtained with the as- 
sumptions made to derive equation H2\ and equation H4\ . In this 
case the dependency of a (or r e ) on the mass M is undetermined (s 
and p cannot be computed). 

However, the second solution yields different values of the co- 
efficients a and b with the same classical assumptions. In this case, 
the dependency of r e or a on M can be computed. Hence the tilt 
of the observed fundamental plane with respect to the virial plane 
could simply be explained in this way whithout relaxing any of the 
assumptions made above. Taking typical ob served values a — 1,1 
and b — 0.3 (e.g. iFraix-Burnet et alj [2010), we find s 2.5 and 
p =i -4. This yields strong dependences of the variables on mass. 
In particular, the r e vs M relation i s steeper than that foun d for in- 
stance in numerical simulations bv lRobertson et ai]( l2006l) . 

As a conclusion, the two possible solutions implied by equa- 
tion d 1 6b , leading to s and p either undefined or relatively high, can 
be debated, but it appears clearly that the assumptions made to de- 
rive the virial plane are a very restrictive case of more general con- 
ditions (equation Jilt ) that lead to a fundamental plane correlation. 
Interestingly, one usual modification to explain the apparent tilt of 
the observed fundamental plane with respect to the virial plane is 
to relax the assumption in equation d 14fc with t + 1 and rely on 
the observations for this value. In the next section we show that we 
can go a step further and use the more general conditions in equa- 
tion dl lb with the sole observations as constraints without any a 
priori input from the virial hypothesis. This allows the exploration 
of wider physics: the virial conditions (equations J 1 2b and l!14b). if 
fulfilled, should be derived from the observations, not taken as an a 
priori. 



3.3 Observational constraints 

Henceforth, let us use the relations in equation dl lb to derive the 
constraints on the functions in equation Q from the observations. 
To derive the three exponents s, p and t, we take the values for a 
and b as provided by the observed fundamental plane correlations, 
and also other observed correlations. 

As mentioned above, the observed fundamental plane yields 
typical values of a = 1.1, b m 0.3 and c = -8.5 (and in = 4.45 in 
the R band). However, these are known to depend on the sample, 
the waveleng th, the redshift and mor e generally on the group of 
galaxies (e.g. lFraix-Bumet et al.l2 010). This strongly indicates that 
the functions of X in equation l|7) and X itself are not universal and 
depend on the population of galaxies and their ass embly histories. 
This c an be seen in the numerical simulations by iRobertson et al.l 
(2006) where the r e vs M relations change according to the se- 
quences and nature of the merging events. 

With the observed values for a, b and c, equation ( fTTT l becomes 



p =: l.l.y + (-2.5f + 5p)0.3 
log Aj = 1 . 1 log A 2 + 0. 3 (2.5 log (*nA\ /A 3 ) + 4.45) 

-8.5 

2.2s+1.5f 



p 

A, = 2 10 



2 A»- 5 



(16) 



The above set of relations is independent of the actual meaning 
of X. The physics can here come to our help to consider different 
possibilities for the confounding parameter X, but here we still only 



use observations (of real objects or synthetic ones from numerical 
simulations) to determine the s hape of the functions of X. 

Let us consider the fig. 5 in Hopk ins et al. I d2008h that gives the 
dependence of r e and a on the mass fraction f starburst of the starburst 
which is an indicator of the level of dissipation in passed merger. 
The dissipation has been found to be a key element to explain the 
very existence of the fundamental plane. In simple words, without 
dissipation, the merger product would not follow a plane, that is to 
say that the correlation of equation ((9) would not exist. Hence if 
we assume that the relations given in their fig. 5 are observations 
(of synthetic populations of galaxies in this case), then we derive 
that p =! -1 and s — 1. We then compute the exponent t from 
equation d 1 6b and find t — 2.\. We thus have 



r 1 

J starburst 

f l 

starburst 
starburst 



(17) 



These functions are in agreement with the expected role of 
dissipation: evacuates angular momentum, concentrates the matter 
in the bulge (r e diminishes), increases the surface brightness, and 
increases the velocity dispersion. 

Let us consider another candidate for X: M BH the mass of the 
central black hole. It can be guessed that it could play a key role 
in the evolution of the prop erties of galaxies in t he course of di- 
versification. On the Fig. 4 in lHopkins et alJd2009l) . there is a clear 
dependence of r e and or on M BH . The relations are not exactly lin- 
ear, but let us assume that it is the case to a first approximation. We 
find that p m 0.63 and s a 0.28, implying from equation d 1 6b that 
t a 0.83. Then 



m 



,0.63 
BH 



(18) 



These values c ompa re very favorably to the results obtained 



by iGultekin et alj d2009h (tr oc M° B t) 5 and L oc M°/ H ) and 



iMarconi & Huntl [20031 )71 oc M°f) for elliptical galaxies. We then 
conclude that the fundamental plane is explained by solely assum- 
ing that the black hole mass evolves in parallel with other properties 
of galaxies, without any direct physical modeling of each relation 
in equation d!8b . 

Finally, it is interesting to remark that we can combine the 
two examples above with X = (fstarbmst, M BH ) like in equation l[5} 
and obtain the conditions equation ©. We then obtain exactly a 
combination of equation J 1 7b and equation J 1 8 b : 



J starburst BH 



0.63 
BH 
■0.2K 



f 1 M' 

J slm-burs! BH 



(19) 



fs 



starburst BH 



This solution is given for illustration only because one might 
question the independence between j 'starburst an d M BH (see below). 

We have thus proven that the fundamental plane correlation 
equation l|9) can be easily and plausibly obtained with confound- 
ing parameters like merger dissipation or central black hole mass, 
without any assumption linked to the virial equilibrium. 



4 EVOLUTIONARY CORRELATIONS 

In the course of diversification, many properties of galaxies change, 
and they tend to statistically change in a more or less monotonous 
way. For instance mass and metallicity are both bound to increase 
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with the complexity of a galaxy assembly history, so that they ap- 
pear to be statistically correlated. It seems difficult to avoid the evo- 
lution to act as a confounding factor. 

We thus propose that the main confounding parameter is X = 
T with T an indicator of the level of diversification, being some- 
thing like an evolutionary clock not necessarily easily related to 
time or redshift. 

In the simple case of a single stellar population, the time 
since formation is naturally a good evolutionary clock for some 
parameters like luminosity L, color and metallicity. However this 
is less obvious for r e or cr. Time since formation can probably be 
as good for a homogeneous population of galaxies if they are not 
affected by significant transformation events (interactions, merg- 
ers...). However, galaxies from the same homogeneous population 
(identical properties) do not form at the same epoch in the Uni- 
verse, so that the time since formation is not related to redshift. In 
addition, galaxies are much more complex than single stellar popu- 
lations, so that the "age" of a galaxy unfortunately is only an aver- 
age over the different stellar populations and does not characterize 
its complete evolutionary stage. 

Considering now X = (1+z), we take from ISagliaetalJfioToh 
that p =! -0.5 and s ^ 0.4. Note that these values are obtained at 
fixed mass, which probably biases significantly the derived evo- 
lutions of r e and cr. Anyhow if a tight fundamental plane corre- 
lation exists with the parameters previously used (those given in 
Sagli a et al. I < l2010h . close to ours, do not change the main result 
here), then equation d!6t yields : 

r e oc (l+z)-°- 5 

cr oc (1+z) 04 (20) 
L oc (1+z) - 25 

Our lu minosity ev o lution is clearly weaker than the one es- 
timated in ISaglia et alj feoich : L oc (1+z) 1 . This discrepancy 
could probably be explained by the various hypotheses they have 
to make to try disentangling all effects in such a multivariate and 
evolutionary context. In particular, the evolution of r e and a are 
computed at fixed mass, possibly introducing interdependencies 
of variables through the M/L ratio. More importantly, the cosmo- 
logical clock (1 + z) is not a good evolutionary clock for a mix- 
ture of different populations of galaxies since they do not evolve 
at the same time, at the same space and along the same paths. 
It is well kno wn that the tilt of the fundamental plane depends 
on the sample dD'Onofrio et alj |2008) or on the group considered 
dFfaix-Burnet et alfeoiofc 

Anyhow, diversification cannot be summarized with only one 
simple property (like redshift or mass) because galaxies are too 
complex objects and do not evolve linearly in a unique way. In some 
diagrams, that is for some set of variables, a particular property 
could crudely depict the general trend of diversification. In the case 
of r e , cr and fi e , and to a first approximation, mass could well repre- 
sent a satisfactory driving parameter for the fundamental plane cor- 
relation, but it is certainly not unique as shown in Sect. 13.31 Since 
it is only approximate, some dispersion is expected. 

A lot of observables evolve with diversification, at least sta- 
tistically, so that we should not be surprised by the many scaling 
relations found for galaxies and the difficulty to pinpoint the driv- 
ing parameters and mechanisms. We also better understand why 
several characteristic parameters (mass, luminosity, metallicity...) 
and also the samples themselves have been found to influence the 
shape of the fundamental plane without providing a clearer picture 
of its origin. 

This might also explain some of the observed dispersion is 



most scatter plots. For instance, it has been found that the disper- 
sion of the fundamental plane st rongly depends on the evolutionary 
group (Fra ix-Burnet et af]|2010h , the correlation equation {9]l even 
not holding in the least diversified groups. Also relations several 
parameters like in equation dl9t are most often wrong because the 
evolutionary clock T makes most variables to be non independent. 
Hence, dispersion may be explained by the statistical (non-causal) 
nature of the correlation and the heterogeneity of the samples as far 
as diversification is concerned. 

Indeed, the evolutionary clock, i.e. the factor X = T, can be 
hidden, not understandable analytically and not directly observable. 
It is more directly related to an evolutionary classification, and is 
a well-known pro blem of comparative methods in phylogeny (e.g. 
lFelsensteirJll985h . 



5 CONCLUSION 

The fundamental plane correlation and similarly scaling relations 
for galaxies can be formalized as confounding correlations. The 
confounding factor X is very probably related to the level of diver- 
sification and may be identified in some cases with some variables 
that trace the global evolution of galaxies. The dependence of the 
variables involved in the observed correlations on X has been as- 
sumed here to be power law functions for sake of simplicity but 
they could be more complicated without changing the result of 
the present paper. In particular, these functions can be multivari- 
ate, with several confounding factors. 

The physics thus should not be invoked to explain causally 
each of the observed correlations, but rather to explain the con- 
founding or evolutionary nature of these correlations. Since the 
galaxy assembly history and transformation processes are complex, 
it is quite improbable that a simple physical theory can yield X 
and the dependence of observables on this parameter. Indeed, these 
functions more probably come from the statistics in the course of 
galaxy diversification. In addition, X is probably hidden, not under- 
standable analytically and not directly observable. It might also be 
different depending on the set of variables and the group of galaxies 
considered. 

Gaining insights on the confounding parameter X requires 
several complementary statistical approaches. A first one is to com- 
bine several scaling relations and several galaxy samples in order to 
identify some common variables that could play the role of the con- 
founding parameter. A second approach is to use numerical simu- 
lations to produce synthetic populations of galaxi es with as many 
assembly configurations a s poss ible (like in e.g. iRobertson et al.l 
|2004 iHopkins et alj|2008l. 120091) . The great advantage here is to 
play with unobservable parameters. Finally, the third approach is 
to group galaxies according to their assembly history. Since the 
confounding parameter X is probably mainly linked to the level 
of diversification, the scaling relations necessarily depend on the 
evolutionary groups. Indeed, the very nature of X and the functions 
characterizing the dependence of the variables on X are expected to 
depend on the assembly history of galaxies. 

Combining evolutionary classifications, numerical simula- 
tions, observed scaling relations and recognizing the latter as evo- 
lutionary correlations, will lead us toward a much better under- 
standing of the history of galaxy formation. Interestingly, all astro- 
physical objects evolve, and evolutionary correlations could prob- 
ably also explain scali ng relations for other ste l lar systems such as 
globu lar clusters (e.g. Misgeld & Hilker 2oTH: iFraix-Burnet et all 
l2009h . 
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